Preprint typeset in JHBP style - HYPER VERSION 



Constraints on single-field inflation with WMAP, SPT 
cn ■ and ACT data — A last-minute stand before Planck 



o 



Cheng Cheng^'^ * Qing-Guo Huang^ t Yin-Zhe Ma 



3,4 t 



^ State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese 
Academy of Sciences, Beijing 100190, China 

^ University of the Chinese Academy of Sciences, Beijing 100190, China 
' Department of Physics and Astronomy, University of British Columbia, Vancouver, 

^ ! V6T IZl, BC Canada 

• I Canadian Institute for Theoretical Astrophysics, 60 St. George Street Toronto, M5S 

Q^' 3H8, Ontario, Canada 

o 

> ' Abstract: We constrain models of single field inflation with the pre-Planck CMB data. 

' The data used here is the 9-year Wilkinson Microwave Anisotropy Probe (WMAP) data, 

South Pole Telescope (SPT) data and Atacama Cosmology Telescope (ACT) data. By 
^ adding in running of spectral index parameter, we find that the is improved by a factor 

. of Ax^ = 4.22, which strongly indicates the preference of this parameter from current 

, data. In addition, we find that the running of spectral index as does not change very much 

T^lj- I even if we switch to different pivot scales, which suggests that the power law expansion 

I of power spectrum is accurate enough till the 1st order term. Furthermore, we find that 

' the joint constraints on r — give very tight constraints on single-field inflation models, 

and the models with power law potential (pP can only survive if 0.9 ^ p ^ 2.1, so a large 
class of inflation models have already been ruled out before Planck data. Finally, we 
. use the fjsfi data to constrain the non-trivial sound speed c^. We find that the current 

^ . constraint is dominated by the power spectrum constraints which have some inconsistency 



with the constraints from /nl- This poses important questions of consistency between 
power spectrum and bispectrum of WMAP data. 
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1. Introduction 

The inflationary model [||, ^, ^ has achieved a great success in modern cosmology, and it has 
been confirmed by many high precision CMB and Large scale structure experiments . 
It provides a good explanation to a series problems such as flatness problem, horizon 
problem, and monopole problem in the standard cosmology scenario. In addition, inflation 
paradigm provides a natural explanation for the origin of primordial perturbations which 
constitute the seeds for the large scale structure we can see today. Therefore, identifying 
the realistic inflation model becomes an important task in observational cosmology. 

Astronomical observations provide a large mount of data to constrain the cosmological 
parameters, especially inflation models. The default cosmology model people always use is 
the "six-parameter" ACDM cosmology model, in which the single-field slow-roll inflation 
(sound speed Cg = 1) is assumed in the model. However, the class of slow-roll inflation 
models already have some weak tension with the observational data. In Ref. it is 
shown that the generic (pP inflation model cannot provide consistent r — Ug values within 
reasonable range of number of e-folds. In addition, WMAP 9-year data Q suggests that 
the local non-Gaussianity has a large positive value, while the orthogonal non-Gaussianity 
is a large negative value, and these values are hardly to be produced in single-field slow-roll 
inflation models. Given these interesting tension between the single-field slow-roll inflation 
model and the current observational data, we would like to explore the possibilities of non- 
trivial sound speed 7^ 1 as well as non-zero running of spectral index dug/dlnk to test 
their consistency with current combination of WMAP 9-year data ACT data |p and 
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SPT data Q. We intend to finish this work right before Planck data release (expected on 
21st March, 2013) in order to make an immediate comparison before and after the Planck 
data. We hope that our work will motivate theorists to explore more phenomena in the 
non-standard inflation model given the tight constraints on single field inflation models. 




Figure 1: WMAP9 temperature data with lensed ACT and SPT data. WMAP9, ACT and SPT 
data are mainly in the range 2 <l < 1000, 300 < 3000, and 700 </ < 3000 respectively. The 
theoretical curve is the lensed CMB power spectrum with WMAP 9-year cosmological parameters 
and the light blue band is the cosmic variance. The Planck data will hirther tighten up the error- 
bars in the middle regime. This figure is reprinted with permission from Mark Helpern. 

This paper is organized as follows. In Sec. we will discuss the model we are focusing 
on, and the data we will use to constrain the models. In Sec. ^, we will present our results 
of fitting. The concluding remarks will be presented in the last section. 



2. Methodology 
2.1 The Model 

We will use standard 6— parameter ACDM model as our basic model. ^. We then allow 
r (tensor-to-scalar ratio), dus/dlnk (running of spectral index), Cg (sound speed for the 
curvature perturbation modes) to be varied since we want to explore the level of constraints 
from these parameters. The sound speed is related to the tilt of tensor power spectrum in 
the general single-field infiation model through Q 



^The free running parameters are {Qth^, ^ch^, Oa, t, Us, As}, which are fractional baryon, cold 
dark matter, and dark energy density, optical depth, spectral index and amplitude of primordial scalar 
perturbation respectively. 
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where the tensor power spectrum is parameterized as 



P,{k) = Mko) ( -] , (2.2) 



(2.3) 



Here kQ is the pivot scale. Thus the tensor to scalar ratio is defined as 

Mko) 

As{ko)' 



Since tensor power spectrum also contribute to CMB angular power spectrum Cf^ on very 
large scales, we will use the CMB temperature angular power spectrum to constrain r and 

Cs. 

In addition, we add the "running of running" parameter which characterizes the run- 
ning of running of spectral index, i.e. 



da, d? 



Thus the scalar power spectrum is parameterized as 



P,{k)=AM^-^J . (2.5) 

Note that once the "running of running" ((3s) is introduced into the model, the running 
of spectral index as becomes a scale-dependent quantity. To remove any ambiguity, we 
need to specify the pivot scale in the power law expansion (Eq. (|2.5D), this is why the Os 
is related to ko. However, if turns out to be less dependent on /cq, it means that the 
truncation till Og is enough (1st order), and there is no need to introduce a higher order 
truncation {/3s)- 

The reason we want to release f3s as the running of running parameter is that SPT data 
m gives a detection of a negative value of the running of spectral index as = dus/dlnk 
at kQ = 0.025Mpc~^. So we would like to add this parameter as a higher order effect to 
monitor any possible "running of running". Even though it has not been detected, it is 
expected to be significantly constrained and is useful for the reconstruction of single-field 
slow-roll inflation 10|. 



2.2 The data 

We will use the most precise class of CMB data up-to-date, which is the combination of 
WMAP 9-year data [§, SPT data || and ACT data g. The temperature angular power 
spectrum from three data sets is shown in Fig. |^. The combined data is named as "CMi? 
data" in the following discussion. We set the maximum /-range of scalar model to be 7000 
i^max = 7000), and maximal tensor /—range to be 3000 {Imax — 3000) in the running of 
MCMC chains. In addition, we add Baryon Acoustic Oscillation data ||Tl[] as well as B.q 



prior from HST project |12| into our data source. In order to explore the variation of 



sound speed, we add constrained /at^ data provided by WMAP 9-year bispectrum into our 
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likelihood. The sound speed Cg is related to the equilateral and orthogonal type of non- 
Gaussianity /nl through Eq.(57) in [Q]. So according to ||7| we assign /j = {f^i^L^ fNL^) ^ 
the data vector, which is 

f^^ = 51 ± 136 (-221 < f^^ < 323 at 95%CL) 
/^'•f = -245 ± 100 (-445 < /^^'^ < 45 at 95%CL). 

Then we use the function (Eq.(58) in|^]) to calculate the best-fit value of Cg, i.e. 

= S fiFijfj - 2 S Fufih + S m,Fr.^Fjjfj (2.6) 



where Fij is the Fisher Matrix 



F 



/ 25.25 1.06 -2.39 

1.06 0.54 0.20 I X 10"^ (2.7) 

V- 



-2.39 2.20 1.00 



For the extended model of ACDM, we will release r and in the CAMB code [13| 
and further modify the code to incorporate running of running parameter (/3s). We run 
CosmoMC Ip^ , |l5| to generate MCMC samples. We will express our results in term of 
best-fit value of marginalized likelihood, as well as la and 2a confidence level (CL) (68.3% 
and 95.4% CL). 



3. Results 

3.1 Canonical single-field slow-roll inflation model (c^ = 1) 
3.1.1 ACDM cosmology model 

We first fix Cs = 1 and investigate the constraints on parameters r and as- The data sets 
we use here are CMB data, BAO and HQ. Here we consider "6— parameter model", 
"6— parameter+r model", "6— parameter -l-Os" model and "6— parameter+r-l-Os" model 
which are expressed as "ACDM", "ACDM-^r", "ACDM+a^" and "ACDM-hr+a^" models 
respectively. 

In Fig. ^, we can see that the likelihood of r shifts a little if we switch Ug on and 
off. The solid line is ACDM-t-r model, and the dotted line is ACDM+r-|-as model. In 
addition, the likelihood becomes broader in ACDM+r-|-as model, and the upper limit is 
also higher. This indicates that without the direct polarization power spectrum, it is hard 
to draw concrete upper limit on the amplitude of tensor mode r, since adding a single 
extra-parameter can greatly broaden the constraint on r. 

Similar thing exists in Fig. ^. The solid line is the ACDM-t-a^, and the dotted one is 
ACDM-|-as-|-r. One can see that the likelihood of is broader if r is released as a free 
parameter. This means that the two parameters have some level of degeneracy, which is 
potentially able to be broken if the future polarization data is added. 

The Fig. |^ shows the likelihoods of n<j for three models. Here we consider all of the three 
models, i.e. ACDM-hr, ACDM-Fa^, KCDM+r+ag. We can see that not only the peak of 
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Figure 2: Likelihood of r in case of as fixed and as as free parameter. 
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Figure 3: Likelihood of in case of r fixed and r free. 

distribution shift, but also the range of confidence level of Ug changes quite a lot in three 
different model: if we add r, the spectral index still prefers a "red" spectrum as < 1, 
but such situation does not exist anymore in the case of ACDM+q^ and ACDM+a^+r. 

Fig. H shows the contours of joint constraints on — Ug in ACDM+a^ model(the blue 
solid curves) and ACDM+r+ag model(the red dashed curves). We can see adding r leads 
to the shift of rig towards bluer region, and the constraints become broader. 

In Table |l], we list the results of fitting by fixing = 1. One can see that by introducing 
as parameter, the really improve significantly (Ax^ = 4.22), indicating that the current 
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Figure 4: Likelihood of ris in different models. 




Figure 5: Joint constraint of Ug and as in ACDM+a^ (blue solid contours) and ACDM+a^+r 
model (red dashed contours). The contours show la and 2a constraints. 



data prefer the inflation model with running of the spectral index. 

In the left panel of Fig. |6|, we compare our joint constraints on r and Ug with the results 
from WMAP 9-year paper Q. WMAP 9-year results used WMAP 9-year data, combined 
with old SPT data, old ACT data, BAO data and Hq data and obtain the black contours 
(Icr and 2a CL). We used the similar combination, except that our SPT and ACT data are 
the corresponding new data sets js], ^. By updating the new data of ACT and SPT, one 
can see that the constraints are tightened up to some extent. This suggests that the new 
SPT and ACT data really provide a large level arm for WMAP9 data, which offer more 
constraining power on small scales CMB angular power spectrum. 

We use our results of joint constraints on plane of r — to discuss its implication for 
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inflation models (Right panel of Fig. 
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Figure 6: Joint constraint on r — rig and comparing with WMAP9+old SPT+old ACT (left panel) 
and model predictions (right panel). Left panel: the red contours are from WMAP 9-year results 
Q which combined WMAFQ+old ACT+old SPT+BAO+i/o, while our constraints are the resuhs 
of WMAPg+new ACT+new SPT+BAO+i7o- Right panel: we consider several typical inflation 
models. (1) chaotic inflation model ||l^ with potential V{(j)) oc (j)^. The solid lines correspond to the 
predictions for different value of p, and the shallow and darker green dashed lines correspond to the 
predictions for = 50 and A = 60 in the models with different power index p. (2) spontaneously 
broken SUSY (SBS) inflation model whose potential is given by V{(t)) = Vo (l + '^^^g) which 
is assumed to be dominated by Vq. (3) mass term (MT) inflation model with potential V{(f>) — 
Vo — ^m^(jp where the mass term is assumed to be subdominant. 

• Chaotic inflation model whose potential is given by V'((/>) cx; (p^. This model 
predicts r = = 1 — where N is the number of e- folds before the end of 
inflation. Given the current constraints on the amplitude of inflation and the "slow- 
roll" parameter, A'^ is around 60 but with some uncertainty of reheating process. Here 
we take the range of 50-60 as the reasonable range of number of e-folds. The region 
between two dashed lines in Fig. ^ indicates the prediction of chaotic inflation. One 
can see that the models with p = 2 and p = 2/3 are disfavored at around 2a 
level, and only the models with p £ [0.9, 1.8] for = 50 or p € [1.5, 2.1] for A^ = 60 
are still consistent with data within 95% CL. 

• Spontaneously broken SUSY (SBS) inflation model with potential V{(p) = Vq fl + cln ^ 
where the potential is assumed to be dominated by Vq and c <C 1. This model pre- 
dicts r = and = 1 — The spectral index in this model is quite large and it is 
disfavored at more than 95% CL. 

• Mass term (MT) inflation model with potential V((/)) = Vq — ^m?(j)'^ where the mass 
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Table 1: Results of fitting by fixing Cs = 1. We set l^^^ax ~ 7000, and Imax — 3000 in the running 
of MCMC chains. 





kCDM 


AC DM + r 


ACDM + as 


ACDM + r + as 


ACDM + r + as+Ps 




0.961 ±0.007 


0.959 ±0.006 


1.018 ±0.027 


1.066 ± 0.040 


1.089 ±0.080 


r{95%CL) 




< 0.12 




< 0.42 


< 0.53 


Ola 






-0.021 ±0.009 


-0.035 ±0.012 


-0.050 ±0.057 


Ps 










0.005 ±0.021 


Best fit -log(Like) 


4921.52 


4921.15 


4917.30 


4916.91 


4917.54 



term is assumed to be subdominant. The tensor-to-scalar ratio and spectral index in 
this model are respectively given by r = and rig = 1 -|- 2?7 where rj = —m?Mp/VQ. 
This model can fit the data very well \{ r] = —0.02. 

3.1.2 Comparison of different pivot scale 

In the former sections, all the fittings are done at pivot scale = 0.002 Mpc~^ and the 
running of spectral index is preferred at more than 2a level. In this section, we investigate 
the distributions of Ug at different pivot scales. We use the model ACDM-|-r-|-as- The solid 
line \sk = 0.002 Mpc"\ and the dotted hne is A: = 0.025 Mpc"^ in Fig. |. It shows 
that when the pivot scale change, the distribution of almost does not change at all. 
This means that the constraints on is not sensitive to the pivot scale you choose, which 
indicates that the truncation of power index expansion (Eq. ( |2.5D ) is accurate enough till 
1st order. 
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Figure 7: The marginalized distribution of running of spectral index as at different pivot scales. 
3.1.3 Running of running of spectral index {(ig) 

We considered the higher order power effect of the primordial power spectrum, and we 
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introduce a new parameter fis to characterize the "running of running" (Eqs. ( |2.4| ) and 
( |2.5| )). Left panel of Fig. ^ shows the joint constraint on and /3s, and the right panel 
shows the marginalized distribution of (ig with a flat prior. We can see that the peak of /3s 
slightly deviates from 0, but is perfectly consistent with zero within la CL. This means that 
the current data do not support the "running of running of spectral index" , and therefore 
the power law expansion of the scalar power spectrum (Eq. (^^)) is accurate enough till 
the as term. This is consistent with what we find in Sec. |3.1.2| . The fitting results are 
shown in Table |I[ 





-0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 



-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 



Figure 8: Left: Joint constraints on as — Ps- Right: Marginalized distribution of fig with Icr CL 
0.005 ±0.021. 

3.2 General single-field infiation Model (c^ free) 




Figure 9: Likelihood of Cs in different datasets. 
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In this section, we release Cs as a free parameter which is constrained by the f^i data from 
WMAP 9-year results [0] ^. Table. ^ shows the best fit (-log(Like)) and confidence level of 
Hs, r and Os of two models. It can be seen that adding Og significantly reduces the best fit 
-log(Like), but enlarge both the confidence interval of r and Ug. 

Fig. ^ shows the distribution of Cg. The dotted line is the likelihood from /tvl data 
while the solid line is marginalized probabilities from the CMB+BAO+i/o+/AfL data. We 
can see that the Jnl prefers a very low value of Cg which is NOT consistent with the CMB 
power spectrum constraints. In addition, the combined constraints are dominated by CMB 
power spectrum simply because the number of CMB power spectrum data is far greater 
than the /nl data. The tension between the /atl data and the CMB power spectrum data 
may have a variety of indications: 

1) it may indicate that the power spectrum and bispectrum data are not consistent with 
each other, which suggests that there are some uncleaned systematics in the data sets; 

2) it may also indicate that the underlying model, i.e. single-field infiation cannot work at 
all when we confront it with CMB power spectrum and bispectrum data. 

In any case, we need to develop a method which can direct relate Cg with power spectrum 
and bispectrum (not just f^i data) and globally fit this parameter by using full spectrum 
of CMB data. Such work is in progress. 

Table 2: Results of fitting with Cg as a free parameter. 





AC DM + Cs 


AC DM + Cs + as 




0.958 ± 0.006 


1.064 ±0.040 


r(95%CL) 


< 0.11 


< 0.40 


as 




-0.034 ±0.012 


Best fit -log(Like) 


4922.14 


4918.17 



4. Conclusion 

In this paper, we combine the most recent pre-Planck CMB data to constrain the infiation 
model parameters. Our data consists of WMAP 9-year data Q, ACT data Q, SPT data 



|5[, Baryon Acoustic Oscillation data [11| as well as Hq prior |12]. We mainly find four 



interesting results from our numerical fitting: 

• if we add in the running of spectral index ag = dug/dln k, the value reduces a lot, 
which indicates that it improves the fit to data very much. 

• By adding in a '3rd-order' parameter, i.e. running of running /3g, we find that current 
data do not support non-zero detection of /3g. In addition, by switching to different 
pivot scales, the constraints on do not vary a lot. These two tests strongly suggest 
that the expansion of the power spectrum is accurate enough till the 1st order (a^ 
term), and there is no observational hint for the higher order scale-dependent terms. 



^ A similar constraint from WMAP 7yr data and Background Imaging of Cosmic Extragalactic Polariza- 
tion (BICEP) experiment was given in [^. 
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• Due to the new ACT and SPT data we used, our constraints on r — is tighter 
than the WMAP 9-year results [^. Our constraints is aheady able to rule out a 
large class of single-field inflation model even before Planck data. We show that the 
single field inflation with power law (jf can only survive if p is in between 0.9 and 
2.1, and Spontaneously broken SUSY (SBS) inflation is ruled out firmly by current 
observational data. 

• We release sound speed Cg as a free parameter, and find that the constraints on Cg from 
fNL data and CMB power spectrum are not consistent with each other. This strongly 
indicates that either there is some unaccounted systematics in the bispectrum data 
that may incur extra-error in the fjsiL estimation, or the model of varying Cg cannot 
work at all given these two datasets. In any case, this motivates us to explore a set 
of formulism that directly compute power spectrum and bispectrum given a Cg value. 

In conclusion, we find that pre-Planck data have already been able to set tight con- 
straints on single field inflation model. But current observational data still leave many 
open questions to be solved. We hope such issues will be resolved when the Planck data 
becomes available in a few days. 
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